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Abstract. Initial chemical composition of stars is, be- 
sides their mass, another key factor in stellar evolution. 
Through stellar lifetimes and impact on radiation output 
and nucleosynthesis of stars it is controlling both the pace 
of evolution of galactic matter/light and changes in their 
integrated observables and spectra. However, because of 
long-time scarcity of homogeneous stellar evolutionary 
data for compositions other than solar, or "normal", all 
the estimates of physical situation and global parameters 
for galaxy collective (radiation density, "global" star for- 
mation rate (SFR) etc.) have effectively remained being 
restricted to the assumption of normal composition and 
the composition-dependence of galactic properties ignored 
even when it should be obviously rather significant, e.g. in 
discussions of the main galaxy-formation epoch between 

z~ l-e-10. 

With our evolutionary synthesis codes we have per- 
formed extensive computations of temporal behaviour of 
photometric and spectral properties of stellar populations 
having different initial Z and varying IMF-s and SFR-s. 
Also, the production of most common nucleosynthesis el- 
ements He, C and O was followed. These computations 
have been performed on the basis of two available but 
different homogeneous multicomposition stellar evolution 
tracks grids by Geneva and Padova groups and the Kurucz 
model atmospheres (Kurucz 1993). 

In given paper next to the discussion of overall effects 
of evolutionary differences we also present and comment 
the normalized per stellar mass unit standard tables of 
the detailed radiation energy output from stellar popula- 
tions with different chemical composition, integrated over 
the whole lifetime of their stars, likewise the tables on 
He, C and O production. The energy production at dif- 
ferent wavelengths per production of unit amount of 16 
and heavy elements are tabulated in the function of the 
"metallicity" Z, the slope parameter n = x+1 (sic!) of the 
initial mass function (IMF) and the lowest stellar initial 
mass boundary, M cut , for stars evolving finally into black 
holes. 
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1. Introduction 

Apart from the initial nucleosynthesis in the first ~ 10 3 
seconds after the Big Bang, what settled the primordial 
abundances of light elements (cf. Wagoner et al. 1967; 
Walker et al. 1991: Smith et al. 1993), the only source 
of the chemical enrichment of the matter in the post- 
recombination Universe has been the nucleosynthesis by 
stars, with its products mixed into the ambient medium 
by stellar ejecta during rapid mass loss stages and through 
stellar winds. This way stellar nucleosynthesis has been 
the major driving force behind all the subsequent evolu- 
tion of stellar systems and the main generator of radiation 
energy and the cosmic background from the formation of 
very first stars up to the present epoch. 

As known, rapid nucleosynthesis in hot matter follow- 
ing the Big Bang settles the initial He content, Yq, and 
produces negligible amounts of 3 D, 6 Li, 7 Li, Be etc. Cos- 
mology cannot tell us about the exact time and rcdshift, 
when the first stars might have appeared in the Universe 
and therefore the first stellar systems started to form, but 
it is believed to have happened at redshifts around z ~ 10 
or even somewhat higher. Originally extremely deficient in 
heavy elements, the chemical composition of the matter of 
these primeval galaxies will evolve rather unsignificantly 
in low-mass low-density objects, but in the high-mass ob- 
jects the enrichment proceeds rapidly. Consequently, the 
chemical composition differences should be of great im- 
portance, when the radiation of young galaxies is being 
analyzed or used to address the cosmological problems. 
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2. The era of first galaxies 

Longstanding observational search for the primeval galax- 
ies, going historically back to a paper by Partridge and 
Peebles (1967), has by now pushed the redshift limit for 
the directly observed starforming objects the farther and 
farther, with the most distant objects presently known 
having z ~ 7. Although seemingly quite rare objects at so 
high rcdshifts, the discovery of quasars at z > 6 (presum- 
ably surrounded by stellar component and fainter galax- 
ies) testifies, that galaxy formation process and birth of 
first stellar populations is in some cases going very inten- 
sively on even in that young the Universe, with age slightly 
less than 1 Gy. For the most distant quasar presently 
known, Jl 148+5251 with redshift of z = 6.419 (Fan et 
al. 2003), Bertoldi et al (2003) give on the basis of mea- 
sured dust emission for its parent protogalaxy an estimate 
of very impressive star formation, with the rate ~3000 
M /yr. Many far-away galaxies have been discovered in 
quasar fields, some of them seem even to be rather evolved. 
For example, deep imaging with Keck telescope of the 
field around PC 1247+3406 (Soifer et al. 1994), having 
z = 4.897, and the z = 3.8 radio galaxy 4C 41.17 (Gra- 
ham et al. 1994) have revealed the hints of a population 
of red galaxies at magnitudes K ~ 21 and fainter. Gra- 
ham et al. point out that the closest companions to the 4C 
41.17 are very red with their emission most likely caused 
by starlight, with lower limit to their age deduced from 
initial starburst 0.5 Gyr, and if sharing the same redshift 
with 4C 41.17, should have formed at z ~8. 

Another but indirect evidence of the ongoing 
star/galaxy formation at these and even higher rcdshifts is 
the absence of continuum absorption shortward of Ly a in 
quasar spectra (Gunn-Peterson effect), what testifies that 
the intergalactic medium has been practically completely 
ionized starting from the very early epochs down to the 
present. This is generally attributed to the stellar/galactic 
processes be they supernova-driven winds either from lu- 
minous (Schwarz et al. 1975) or low-mass (Tegmark et al. 
1993) high-z young galaxies, cosmic rays (Ginzburg and 
Ozcrnoy 1965; Nath and Biermann 1993) or the photoion- 
ization by stellar (Songaila et al. 1990; Miralda-Escude 
and Ostriker 1990) or quasar (Arons and McCray 1969) 
UV-radiation. Of course, in reality all these sources get 
combined in producing the observed outcome. More exotic 
explanations have also included cosmological blast waves 
(Ostriker and Ikeuchi 1983), decaying massive neutrinos 
(Sciama 1990) etc. 

First discovered by Gunn and Peterson (1965), the 
constraints to the effect got further developed by a num- 
ber of observers, e.g. by Giallongo et al. (1994). They de- 
duced, measuring the average depression between Lyman 
absorption lines of the quasar BR 1202-0725 (z em = 4.695) 
the value of optical depth r < 0.02 ± 0.03 at z ~ 4.3 
and, taking the quasars to be the prime ionizators of 
the gaseous medium and maximizing their known statis- 



tics, the estimate for the density of the IGM Qigm ~ 0.01. 
Juxtaposition of that result with the mean baryonic den- 
sity deduced either from HI Lya absorption at z = 1.9, 
tt b = 0.044 (Tytler et al. 2004), or from D/H measure- 
ments in nearly pristine gas clouds on the basis of conven- 
tional Big Bang nucleosynthesis, fib — 0.0214+0.0020 h^ 2 Q 
(Kirkman et al. 2003, for the popularly accepted value 
H = 71 kms^Mpc- 1 it gives n b = 0.0425 ± 0.004) in- 
dicates, that even at these very early times the dominant 
part of the baryonic matter had been confined into dense 
substructures, what occupied a rather limited fraction 
of the total volume. Vice versa, identifying these dense 
structures at least partly being young/primeval galaxies 
contributing to the ionizing UV-radiation at these and 
higher rcdshifts, it would be possible to satisfy the balance 
Qigm + Qgal ~\~ &hi = ^6, with £l ga i and £Ihi denoting the 
densities of galactic matter and neutral hydrogen clouds, 
respectively. 

3. Spectrophotometric models 

Integrated spectra are the main sources of information 
on the stellar content and interstellar medium of external 
galaxies farther than the Local Group. Being primarily 
built up by thermal radiation of stellar component, the in- 
tegrated spectrum depends on its age, star formation rate, 
stellar mass spectrum and chemical composition distribu- 
tion, but may also be heavily influenced by the reradiation 
and absorption by dust and heated gas in starbursts. Ow- 
ing to specific frequency behaviour, the nonthermal com- 
ponent can be easily recognized , all the other factors work 
together in build-up of the spectrum in a complex manner. 

Evolutionary modelling of galactic spectra was started 
by Bruzual (1983) more than two decades ago, but re- 
mained for most of that period mainly constrained to so- 
lar composition, due to both lack of homogeneous stellar 
evolutionary track sets for different metallicities,as also 
underlined in the abstract, as well as spectral libraries 
for nonsolar compositions. They became gradually avail- 
able to the midst of 1990's. Presently there do exist two 
consistent track sets with sufficient range of metallicities 
- set for 6 compositions produced by Geneva group (cf. 
Meynet et al. 1994, and Charbonncl et al. 1996, with ear- 
lier references in these; a sixth, high-metallicity Z = 0.10 
composition got added to the set later by Mowlavi et al. 
1998) and Padova tracks ( cf. Girardi et al. 1996, and 
references therein, Girardi et al. 2000 for stars of lower 
masses), the latter covering with 8 compositions practi- 
cally all the observable and theoretically useful metallic- 
ity interval from Z = 0.0001 until Z = 0.10. Both yet still 
pose some smaller problems for applications(ors) - Geneva 
tracks do not include for 4 non-solar compositions (except 
Z = 0.001) past-RGB stages of evolution of low-mass stars 
M < 1.7 M Q , so these have to be added from other (and 
inhomogeneous) sources, it also has shorter range towards 
extremely metal-deficient abundances. Plus, the location 
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of stars with masses M < 0.8 M0O11 and subsequent evo- 
lution from the zero-age main sequence is not available, it 
has to be taken elsewhere. Padova tracks are relatively free 
of so basic problematics, except of lack of stars M > 9 M Q 
in Z = 0.10 subset, some minor irregularities arise also for 
Z = 0.001 composition because they were computed with 
different radiative opacities. 

The only set of model atmospheres, suitable for evo- 
lutionary spectral syntesis, has been produced by Kurucz 
(1993), covering spectral range 90 A- 160 ^m for twenty 
compositions. It is known that since Kurucz models do 
not contain bands of some molecules, in the cool star re- 
gion T c fr< 4500°K improvements are needed, but being 
developed by a group associated with Uppsala, such more 
sound atmospheres for cool stars have not became avail- 
able yet. 

At the end of 1970-s and beginning of 1980-s, I de- 
veloped a code package for computions of photometric 
evolution for stellar systems, using theoretical tracks and 
observational color indices and bolometric corrections for 
stars of different temperatures and luminosity classes. 
This package was somewhat complemented at the begin- 
ning of 1990-s as to allow on the basis of stellar atmo- 
spheres to compute spectral distributions for single popu- 
lations as well as for the case of continous star formation in 
model galaxies and to compute the photometry for filters 
within their passbands. With this package, I have per- 
formed a wide set of computations of photomctrical and 
spectral population models with different chemistry and 
input physics (inch also variations of stellar track sets). 
We have used both track sets in these studies of spectral 
evolution of stellar populations in function of their initial 
metallicity and star formation parameters. No gas emis- 
sion/dust absorption has been considered in the models, 
since it is highly individual for galaxies. 

Unfortunately, being made about ten years ago, the 
results of these computions have not been published in 
due fashion. On the basis of Padova tracks and isochrones 
a standardized grid of spectral data files has been released 
(Traat 1996) on a CD-ROM, produced as a collective effort 
to compile a database for galaxy evolution modeling. The 
short description of that grid has been published in the 
paper, accompanying that CD-ROM and describing its 
content (Leitherer et al. 1996). 

That grid includes 1008 models and has the widest 
practically useful range of star formation prescriptions. 
It covers all the Padova set abundance range, i.e. eighth 
mctallicities Z = 0.0001, 0.0004, 0.001, 0.004, 0.008, 0.02, 
0.05 and 0.10. For each Z value that dataset includes 

— six power law IMF-s of different slopes 1.6 -f- 3.5, 

— corresponding single-generation populations (instanta- 
neous formation, " initial starbursts" ) 

— populations with continuous/continuing star forma- 
tion for those IMF-s with 4 SFR index values s = 
(constant SFR), s = 1 (exponentially declining), s = 
1.5 and 2 (initially faster, later slower than the expo- 



nential SFR) and 5 star-formation timescales 0.2, 1, 2, 
5 and 15 Gy. 

Spectra were presented for 50 ages in the case of starburst 
populations and 20 ages for populations with continuous 
star formation. Ages range from 0.003 Gy to 20 Gy. 

SFR has been parametrized by a power of the gas 
volume density (as introduced by Schmidt (1959)), with 
index s and time-scale to. In this context, the single- 
generation ("initial-burst") populations with their inde- 
pendency on the SFR/its power index form the limiting 
"starburst" t — case. 

4. Discussion of the metallicity-dependent spectral 
data 

In the context of model stellar populations and their 
comparison with observations, they fall broadly into two 
cathegories, the passive (initial starburst) and active star- 
forming populations, the first formed in simultaneous 
burst and lacking any subsequent star formation, the sec- 
ond characterized by ongoing star formation. Later evo- 
lution of passive populations is just pure stellar physics 
except the number distribution of stars by mass (IMF), so 
under assumption of perfectness of our knowledge of stel- 
lar evolution it depends solely on chemical composition 
and IMF parameters. Active populations, being a super- 
position of infinite number of passive ones, additionally 
need for their characterization the SFR with its parame- 
ters. 

Spectral evolution of passive, initial-burst populations 
computed with six compositions available in the Padova 
stellar tracks grid is represented on Fig-s 1-2. The method 
used for their derivation consists of the integration of their 
output over stellar mass along the continuous isochrones, 
also, to ensure a perfect "coherence" in metallicity Ku- 
rucz model atmosphere flux tables were interpolated to 
the values of actual track set metallicities before assigning 
spectra to the stars. 

General pattern of the behaviour of the passive popu- 
lations is as follows. Were we having included a pre-main- 
sequence stage of stellar evolution with time-count set-in 
at the Hayashi border, with the result of massive stars 
reaching the main sequence shortly after formation and 
less massive ones only in 0.1-^1 Gy, we would have started 
the population with lower initial blue luminosities, with 
their maximal value reached the later the greater the IMF 
slope, n. As it is not the case here, the populations con- 
sidered have a maximum luminosity at all wavelengths at 
the onset, its absolute value dropping with the lowering of 
the number of massive stars in the IMF(= increase of n) 
and the monotonous decline of luminosity with time. The 
influence of n on the spectrum is the most easily under- 
standable factor, since its growth causes the mechanical 
decline of the massive end of the IMF and, respectively, 
the overall mechanical decrease of the spectrum level to- 
wards shorter wavelengths and rise towards infrared, with 
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Fig. 1. Evolution of spectra of initial burst stellar populations with different metallicities. The Z = 0.0004 spectrum has a 
correct location in log E\, all the others have shifted downwards by additional 0.5 dex in respect of the previous one, with 
maximum shift —2.5 in the Z = 0.05 case. 



the midst of the optical region being the approximate bor- 
derline (Fig. 3b). Both Z (Fig-s 1-2) and the age (Fig. 3a) 
work in the same direction, the short wavelength region 
becoming progressively eroded with their growth, and the 
more, the shorter the wavelength. However, with the rise 
of Z and the enhancement of absorption in metallic lines in 
UV the bulk of energy reradiated in the optical and near- 
infrared is increasing with the definite net result that the 
spectrum level at A > 1 /jm is for all ages progressively 
higher for higher metallicities. So a sufficiently long base- 
line for spectral or multicolor photometry coverage (plus 
some improvement of precision of both theory and obser- 



vations) could probably resolve the age-metallicity degen- 
eracy, I think a theoretical calibration might be attempted 
after the release of improved cool-star atmospheres. 

There are yet some more noteworthy details pointed 
out on Fig-s 1-2 like the extremely rapid time evolution of 
the He° continuum below 504Aand bit shallower change of 
the Lyman jump, making both of them theoretically good 
age estimators for the extremely young stellar populations. 
Fig. 3a is also embedding a reference to a specific pecu- 
liarity of the Padova track system: the high-metallicity 
population includes a hot h&B-manque (=naked stellar 
core) stellar evolution phase, stars in this stage are the 
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Fig. 2. (continuation of Fig. 1). Spectra of the same initial burst stellar populations at very old ages 5 and 15 Gy. Shifts on 
graphs are as indicated in caption to Fig. 1. 



cause of the UV-fiux upturn at ages > 10 Gy. This in- 
teresting result is in full detail covered by the authors of 
tracks (Bressan et al. 1994). 



Spectral evolution of active populations is governed to 
a large extent by the SFR shape (i.e. power index s in 
the case here) and star- formation timescale, to- The red- 
dest limiting case for both the spectrum and colors being 
the passive io=0 population, the active population with 
the given s forms during its youth stars the slower the 
larger to value is, but has higher star-formation intensity 
and bluer colors at later epochs. The limits to what ex- 
tent the selection of the functional form parameter s may 
influence the spectrum are approximately the same, the 
dependence of spectral evolution on s being for one pa- 
rameter set illustrated on Fig. 4a (s=0 =constant, s=l = 
exponentially declining star formation rate). The scale of 
impact depends also on the age of the population, in the 
s=0, t =l Gy model the star formation has had an dura- 
tion of first 1 Gy and not been working for 4 Gy-s until our 
time of observation, it makes the population spectra red 
and close to that of a passive population at the same age, 
5 Gy. Right panel depicts the temporal evolution of an ac- 
tive stellar population with a quite extended timescale, so 
although both its colors and spectrum are slowly evolving 
to red, no impressive changes next to its overall absolute 
decrease are evident, even the UV-branch keeping nearftp 
ly the same level in respect to the optical region, as a 
result of the continuing star formation at the levels not 
drastically differing from those in earlier stages. 



5. Metallicity and radiation output in specific 
wavelength ranges 

Metallicity growth affects the stellar evolution in two 
ways: first, making stars cooler and dimmer, secondly, re- 
distributing their flux through opacity growth at short 
wavelengths towards longer wavelengths. So the summatic 
result of Z rise on the composite spectrum of a stellar pop- 
ulation is the progressive erosion of flux in ultraviolet re- 
gion, and the faster, the shorter the wavelength. However, 
with the growth of Z and the enhancement of absorption 
in metallic lines in UV the bulk of energy reradiated in 
the optical and near-infrared is increasing with the defi- 
nite net result that the spectrum level at A > 1 /xm is for 
all population ages progressively higher for higher metal- 
licities. 

Table 1 gives a general quantitative review of the ex- 
tent of metallicity effects in different spectral regions for 
populations with different mass function slopes n. Time is 
eliminated by integration over the lifetimes of stars, pop- 
ulation mass is scaled to the unit amount of mass in very 
massive stars M > 10 M Q , Geneva set of tracks (5 compo- 
sitions, the Z = 0.10 subset was not included) was used. 
These data testify, that composition-caused flux changes 
can be rather impressive, extending to factor of 10 in far- 
UV and ~ 2 in infrared. 

As to illustrate metallicity effects graphically, we also 
provide a couple of examples on Fig. 5, computed with 
Padova tracks. On the left panel of this figure we have 
plotted computed spectra of coevally formed stellar pop- 
ulations (so-called initial burst populations) with 8 dif- 
ferent initial mctallicities Z, having the same ages 1 Gy. 
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Fig. 3. Left panel: Age dependency of high-metallicity initial-burst population spectrum. Right panel: Impact of the IMF slope 
n on the total radiation energy output of a stellar population, integrated over its lifetime. Given panel is based on the Geneva 
tracks, distributions are normalized at 5550 A. Arrows indicate the direction of growth of n, as to facilitate identification of 
curves. 




Fig. 4. Left panel: Spectral evolution of a stellar population with continuing star-formation, having time-scale to = 5. In 
the given s=2 case half of the matter will during to be converted to stars. Right panel: Influence of the SFR time behaviour 
differencies, as introduced by its power index s values, on the spectrum of a Z = 0.05 population with continuous star formation. 



Peeter Traat: Evolution of stellar populations with different Z and chem. enrichment 



7 



Table 1. Time-integrated energy fluxes of coeval stellar populations with different IMF slopes and metallicities, per 1 
M Q in massive stars 10 -f- 120 M . Flux units 10 52 erg, flux fractions f ; computed for intervals 90.9-395 A, 395-995 A, 995-1995 
A, 1995-3590 A, 3590-7490 A, 0.749-4.99 fim and 4.99-160 /im. 
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Fig. 5. Composition dependence of spectra of stellar populations: coevally formed young stellar populations (age 1 Gy, left 
panel) and old populations with constant SFR (age 10 Gy, right panel). 



Such comparatively blue spectra are typical for medium- 
aged star clusters of assumed chemical compositions, or 
might be relevant to different parts inside young ellipti- 
cal galaxies with ages somewhat exceeding 1 Gy. The flux 
of models is scaled to the unit mass in luminous stars 
with masses M > 0.6 M©, the IMF in these graphs is a 
power-law with slope n = —2.35 ("Salpeter" value). The 
most metal-deficient, Z = 0.0001 spectrum has a correct 
location in log E\, all the others have been successively 
shifted downwards by additional 0.5 dex, with maximum 
total shift —3.5 for the Z = 0.10 case. On the right panel 
the spectra of old, 10 Gy stellar populations, are plotted, 
in which star formation is continuous and proceeds with 
a constant, time- independent intensity over the eon to. 
Given case might be considered as a kind of approxima- 
tion to late spirals or irregulars, since in many of these 
actual SFR-s do not seem to significantly differ from the 
mean average over their past. Notice, however, that due 
to the actively continuing star formation the ultraviolet 
flux of models keeps sizable. The growth of metal content 
sharply reduces the flux at shorter wavelengths, as also on 
the left-panel plot. 

6. Nucleosynthesis 

The tabulated nucleosynthesis of He, C and O will be 
available in replacement paper to this first draft! 
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